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Abstract 

Large entropy fluctuations in an equilibrium steady state of classical me- 
chanics were studied in extensive numerical experiments on a simple 2-freedom 
strongly chaotic Hamiltonian model described by the modified Arnold cat 
map. The rise and fall of a large separated fluctuation was shown to be de- 
scribed by the (regular and stable) "macroscopic" kinetics both fast (ballistic) 
and slow (diffusive) . We abandoned a vague problem of " appropriate" initial 
conditions by observing (in a long run) spontaneous birth and death of arbi- 
trarily big fluctuations for any initial state of our dynamical model. Statistics 
of the infinite chain of fluctuations, reminiscent to the Poincare recurrences, 
was shown to be Poissonian. A simple empirical relation for the mean period 
between the fluctuations (Poincare "cycle") has been found and confirmed in 
numerical experiments. A new representation of the entropy via the variance 
of only a few trajectories ("particles") is proposed which greatly facilitates 
the computation, being at the same time fairly accurate for big fluctuations. 
The relation of our results to a long standing debates over statistical "irre- 
versibility" and the "time arrow" is briefly discussed too. 
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1 Introduction: 

macroscopic vs. microscopic fluctuations 

Fluctuations are inseparable part of the statistical laws. This is well known 
since Boltzmann. What is apparently less known are the peculiar properties of rare 
big fluctuations (BF) as different from, and even opposite in a sense to, those of 
small stationary fluctuations. In this paper we consider the simplest type of chaotic 
dynamical systems, namely a finite-freedom Hamiltonian system which admits the 
(stable) statistical equilibrium (SE). This class of dynamical models is still popular 
(since Boltzmann!) in debates over the dynamical foundations of statistical mechan- 
ics (see, e.g., "Round Table on Irreversibility" in 0, and 0). 

A fairly simple picture of BF in such systems is well understood by now, though 
not yet well known. To Boltzmann such a picture was the basis of his fluctuation 
hypothesis for our Universe. Again, as is well understood by now such a hypothesis 
is completely incompatible with the present structure of the Universe as it would 
immediately imply the notorious "heat death" (see, e.g., |§). For this reason, one 
may even term such systems the heat death models. Nevertheless, they can be and 
actually are widely used in the description and study of local statistical processes in 
thermodynamically closed systems. The latter term means the absence of any heat 
exchange with the environment. Notice, however, that under conditions of the expo- 
nential instability of motion, typical for chaotic systems, the only dynamically closed 
system would be the "whole Universe". Particularly, this excludes the hypothetical 
"velocity reversal" also popular in debates over "irreversibility" since Loschmidt (for 
discussion see, e.g., j4|). 

In any case, dynamical models with SE do not tell us the whole story of either 
the Universe or even a typical macroscopic process therein. The principal solution 
of this problem, unknown to Boltzmann, is quite clear by now, namely, the "equilib- 
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riumfree" models are wanted. Various classes of such models are intensively studied 
today. Moreover, the celebrated cosmic microwave background tells us that our 
Universe was born already in the state of a heat death which, however, fortunately 
to us all became unstable due to the well-known Jeans gravitational instability || . 
This resulted in developing of a rich variety of collective processes, or synergetics, 
the term recently introduced or, better to say, put in use by Haken |5[]. The most 
important peculiarity of such a collective instability is in that the total overall re- 
laxation (to somewhere ?) with ever increasing total entropy is accompanied by an 
also increasing phase space inhomogeneity of the system, particularly in tempera- 
ture. In other words, the whole system as well as its local parts become more and 
more nonequilibrium to the extent of the birth of a secondary dynamics which may 
be, and is sometime, as perfect as, for example, the celestial mechanics (for general 
discussion see, e.g., @, |7|) Hi)- 

We stress that all these inhomogeneous nonequilibrium structures are not BF 
like in SE but are a result of regular collective instability, so that they are imme- 
diately formed under a certain condition. Besides, they are typically dissipative 
structures in Prigogine's term |J due to exchange of energy and entropy with the 
infinite environment. The latter is the most important feature of such processes, 
and at the same time the main difficulty in studying the dynamics of those models 
both theoretically and in numerical experiments which are so much simpler for SE 
systems. 

In the latter BF consists of the two symmetric parts: the rise of a fluc- 

tuation followed by its return, or relaxation, back to SE (see Figs.l and 2 below). 
Both parts are described by the same kinetic (e.g., diffusion) equation, the only 
difference being in the sign of time. This relates the time-symmetric dynamical 
equations to the time-antisymmetric kinetic (but not statistical!) equations. The 
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principal difference between the both, some times overlooked, is in that the kinetic 
equations are widely understood as describing the relaxation only, i.e. increase of 
the entropy in a closed system, whereas in fact they do so (at least in SE) for the rise 
of BF as well, i.e. for the entropy decrease. All this was qualitatively known already 



to Boltzmann |10| . The first simple example of symmetric BF was considered by 
Schrodinger |TT|. Rigorous mathematical theorem for the diffusive (slow) kinetics 
was proved by Kolmogorov in 1937 in the paper entitled "Zur Umkehrbarkeit der 
statistischen Naturgesetze" ("Concerning reversibility of statistical laws in nature") 



| I~2|j (see also [13|]). Regrettably, the principal Kolmogorov theorem still remains 
unknown to both the participants of hot debates over "irreversibility" as well as the 
physicists actually studying such BF (see, e.g., |T^|). 

By now, there exists the well developed ergodic theory of dynamical systems (see, 
e.g., |15|]). Particularly, it proves that the relaxation (correlation decay, or mixing) 
proceeds eventually in both directions of time for almost any initial conditions of 
a chaotic dynamical system. However, the relaxation does not need to be always 
monotonic which simply means a BF on the way, depending on the initial conditions. 
To get rid of such an apparently confusing (to many) "freedom" we take a different 
approach to the problem: instead discussing the "true" initial conditions and/or a 
"necessary" restriction of those we start our numerical experiments at arbitrary ini- 
tial conditions (most likely corresponding to SE), and do observe what the dynamics 
and statistics of BF are like. Of course, such an approach is based on a fundamen- 
tal hypothesis that all the statistical laws are contained in, and can be principally 
derived from, the underlying fundamental (Hamiltonian) dynamics. To the best of 
our knowledge, there is as yet no contradiction to this principal hypothesis. Notice, 
however, that such an approach can be directly applied to the fluctuations in finite 
systems with statistical equilibrium only (for discussion see [Q and ]16[). In such, 
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and only such, systems infinitely many BF grow up spontaneously independent of 
the initial conditions of the motion. This is similar to the well-known Poincare 
recurrences (for farther discussion see Section 4 below). 

In spite of the essential restrictions the simple SE models allow us to better 
understand the mechanism and role of BF in the statistical physics. Besides removal 
of a vague problem of the initial conditions such models help a lot in clarifying 
the relation between macroscopic and microscopic description of chaotic systems. 
Particularly, spontaneous rise of a BF out of SE is a macroscopic event as well as 
its subsequent relaxation back to SE, even in a few-freedom system. Like other 
macroscopic processes the BF are not only perfectly regular by themselves but also 
surprisingly stable against any perturbations, both regular and chaotic. Moreover, 
the perturbations do not need to be small. At first glance, it looks very strange in 
a chaotic, highly unstable, dynamics. The resolution of this apparent paradox is in 
that the dynamical instability of motion does affect the BF instant of time only. As 
to the BF evolution, it is determined by the kinetics whatever its mechanism, from 
purely dynamical one, like in model (2.2) we will use in this paper, to a completely 
noisy (stochastic) one. As a matter of fact, the fundamental Kolmogorov theorem 



12j is related just to the latter case but remains valid in a much more general 



situation. Surprising stability of BF is similar to the less known conception of 



robustness for the Anosov (strongly chaotic) systems [17 whose trajectories get 
only slightly deformed under a small perturbation (for discussion see [[|). 

In the present paper we consider a particular type of BF which is characterized 
by a large concentration of "particles" in a small phase space domain of a dynam- 
ical system. In other words, "our" fluctuations are localized in phase space and 
separated in time. A more accurate definition of these fluctuations will be given 
below in Section 3 (see Eq.(3.6)). The same fluctuations in a stochastic model (with 
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noise) were studied in detail in |L4). There exist, of course, many other fluctuations 
with their own peculiarities (see, e.g., flT3|]). The primary goal of our studies was 
the macroscopic kinetics of big fluctuations on the background of small stationary 
microscopic fluctuations. A brief outline of our results was presented in |TF 

2 A Hamiltonian model: 

most simple but strongly chaotic 

The systems with SE allow for very simple models in both the theoretical analysis 
as well as numerical experiments which the latter are even more important for us. 
In the present paper we will use one of the most simple and popular model specified 
by the so-called Arnold cat map (see ^(J ) : 



p = p + x mod 1 , j\ 

x = x + p mod 1 

which is a linear canonical map on a unit torus. It has no parameters, and is chaotic 
and even ergodic. The rate of the local exponential instability, the Lyapunov expo- 
nent A = In (3/2 + a/5/2) = 0.96, implies a fast (ballistic) kinetics with relaxation 
time t r ~ 1/A ~ 1. Throughout the paper t denotes the time in map's iterations. 
A minor modification of this map: 

p = p + x - 1/2 mod C , . 

x = x + p — C/2 mod 1 

where C is a circumference of the phase space torus allows us to study both a fast 
(exponential) ballistic kinetics ( for C — 1) as well as the slow (diffusive) relaxation 
in p (for C ^> 1) with characteristic time t p ~ C 2 /4D P 3> 1 where D p = 1/12 is the 
diffusion rate in p. In contrast to slow diffusion in p, the relaxation time in x does 
not depend on C (t r ~ 1) so that subsequent values of x are practically uncorrelated. 
The map (2.2) has the (unstable) fixed point at x — xq = 1/2 , p = po = C/2. 
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A convenient characteristic of BF size is rms volume (area) in 2D phase space 
(x,p) 

a(t) = a p (t) ■ a x (t) (2.3) 

occupied by a group of N trajectories (particles). In ergodic motion at equilibrium 
a = a = C/12. Due to a severe restriction to small iV<,iO in the numerical 
experiments (see below) we have to use simple (average) characteristics only like 
(2.3). On the other hand, these are just the macroscopic variables we are interested 
in. 

Below we also restrict ourselves to a particular case of BF with the fixed pre- 
scribed position in the phase space: 

1 C 

Xfi = x = - , p f i = p = — (2.4) 

Then, the variance of phase space size v = a 2 = a 2 ■ a 2 is determined by the relation 

°l = (p 2 ) - Pi o* = (* 2 > - 4 (2.5) 

where brackets (...) denote the averaging over N trajectories. In ergodic motion 
at equilibrium v = vse = C 2 /12 2 . In what follows we will use the dimensionless 
measure v = v/vse v, and omit tilde. In diffusive approximation of the kinetic 
equation the variable v(t) is especially convenient as it is varying in proportion to 
time. Moreover, v — > v p in this case because of a quick relaxation v x — > 1 in x. 

In all advantages of variable v its relation to the fundamental conception of 
entropy is highly desirable. The standard definition of the entropy, which can be 
traced back to Boltzmann, reads: 

S = -(\nf(x,p)) +S (2.6) 

where f(x,p) is a coarse-grained distribution function, or the phase-space density, 
and 5o an arbitrary constant to be fixed later. Notice that the distribution calculated 
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from any finite number of trajectories is always a coarse-grained one. However, the 
direct application of Eq.(2.6) requires too many trajectories, especially for BF of a 
small size. Nevertheless, just in the latter case, which is the main problem under 
consideration, we have found a simple approximate relation 

S(t) « ±\nv(t) (2.7) 



which, at least, gives a rough estimate for the entropy evolution fllB] . Moreover, if 
the distribution is Gaussian 

f(x,p) -> f(p) = = (2.8) 

V Zttv 

the estimate (2.7) becomes exact as it is directly followed from the definition of the 
entropy (2.6). Comparison of both relations for the entropy will be considered by 
the end of next Section for a typical BF. 

A great advantage of (2.7) is in that the computation of S does not require very 
many trajectories as does the distribution function. In fact, even a single trajectory 
is sufficient as Fig.l in |16| and below demonstrate ! 



A finite number of trajectories used for calculating the variance v is a sort of the 
coarse-grained distribution, as required in relation (2.6), but with a free bin size 
which can be arbitrarily small. 

Now we may turn to the numerical experiments. 

3 Macroscopic kinetics: 

complete, regular, and stable 

In this Section we consider the regular BF kinetics. The data were obtained from 
simultaneous running of N trajectories for very long time in order to collect suf- 
ficiently many BF for reliable separation of the regular part of BF, or the kinetic 



subdynamics in Balescu's term (see [21j and references therein), from the stationary 



fluctuations. The separation was done by the plain averaging of individual t>j values 
(i = 1, n) over all n BF collected in a run. 

The size of BF chosen for further analysis was fixed by the condition that current 

v(t) < v b (3.1) 

at some time instant t « tj, the moment of BF. Here a prescribed is the main 
input parameter of the run. This condition determines, in fact, the border of the 
whole fluctuation domain (FD): < v < Vj,- 

The event of entering FD is the macroscopic "cause" of BF whose obvious "effect" 
will be subsequent relaxation to the equilibrium. However, and this is the main point 
of our study, the second "effect" of the same "cause" was preceding rise of BF in 



apparent contradiction with the "causality principle" (for discussion see |16j and 
Section 4 below). In any event, the second effect requires the permanent memory of 
trajectory segments within some time window w which is another important input 
parameter of the run. 

The exact procedure of data processing during the run was as follows. Starting 
from arbitrary (random) initial conditions the selection rule (3.1) is checked at each 
iteration. Suppose, it is satisfied at some instance ij n when the bundle of trajec- 
tories enters FD. In the first approximation we could consider it as the fluctuation 
maximum (or variance minimum): tj = t in where sub i is the number of current 
fluctuation in a run. However, such a simple procedure would cause asymmetry 
with respect to t — ti. A better choice were the rule: ij = (t in + t out )/2 where t out is 
the exit moment of time from FD. Instead, we have accepted the following, a more 
complicated procedure which, as we hope, better restores the true BF symmetry. 
Starting from the moment ti n we search for the minimum of v(t) inside a rather large 
interval ti n < t < ti n + w. If a minimum is found at some t = t min , we check that it 
is the minimum inside the next interval t m i n < t < t m i n + w too. If so, we identify 



this minimum with the BF top, and set: t { = t min , otherwise we put t min to the 
time of better minimum and repeat the last step again. Obviuosly, the parameter 
w should be small compared to (P), the mean period of BF, but sufficiently long 
for the trajectory to leave the FD (3.1). Typically, w ^C 2 , the total diffusion time, 
was chosen. After fixing the current U value the computation within the interval 
U < t < ti + w had been completed, and only then the search for the next BF was 
continued. 

As was already mentioned above, there are two rather simple limiting cases of 
generally very complicated kinetics, namely the fast (ballistic) and slow (diffusive) 
ones. An example of both in one run for N = 1 (!) is presented in Fig.l for two 
fluctuations of different size. In this case the general condition (3.1) was checked 
separately for p and x: 

v p (t) < v pb and v x (t) < v xb (3.2) 

with v pb = v xb ~ 1(T 5 and v b = v pb ■ v xb ~ 1(T 10 . 

The fast part of kinetics is approximately described as 

v(t) « v(0) ■ exp (4At) (3.3) 

where r = t—U, A the Lyapunov exponent (Section 2), and where v(0) ~ 10~ 13 is the 
minimal variance averaged over all n fluctuations observed in the run. Notice that 
the latter value is considerably smaller than the border v b ~ 10 -10 . This is because 
of penetration of trajectories into FD. Interestingly, the ratio v b /v(0) = 2000 is the 
same for both runs in Fig.l 

Surprisingly sharp crossover to diffusive kinetics, clearly seen in Fig.l, is related 
to the dynamical scale of diffusion which corresponds to a certain size Vd of the 
increasing variance at which the exponential growth stops. Roughly, it occurs at 
time r = Td when already \x — xq\ ~ \p — Po\ ~ 1/2, whence v x d ~ 12/4 = 3 and 
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v P d ~ 3/C 2 . Hence, we can characterize the dynamical scale as 

v{r d ) = v d = F d ■ v pd ■ v xd = — , r d = — (3.4) 

where F d is an empirical factor, and r d is found from Eq.(3.3). From data in Fig.l 
the dynamical scale v d 0.015 independent of v b which gives the empirical factor 
F d » 1/3. 

In the diffusion region (v > v d ) the initial kinetics is described by a simple 
relation for the free diffusion (see Section 2): 

v(r) » +v d , r d < r « C 2 (3.5) 

which is also shown in Fig.l. It includes two corrections, r d and v d , due to the 
exponential ballistic kinetics. The first one with opposite signs for two symmetric 
parts of the fluctuation takes account for the "lost" time after (or prior to) anti- 
diffusion (diffusion) while the second correction describes a finite fluctuation size at 
the crossover from (to) the diffusion. The mean empirical r d = 7, used in Fig.l, 
is close to the value r d = 6.5 found from Eq.(3.4) with another empirical quantity 
v d = 0.015. 

The large ratio 

B = ^ » 1 (3.6) 

of the mean fluctuation period (P) to a characteristic time of diffusion relaxation 
(see Eq.(3.5)) is the definition of big fluctuation. It guaranties the separation of 
successive fluctuations in time. 

Now we turn to the main subject of our studies, the purely diffusive kinetics of 
big fluctuations. To this end we, first, get rid of re-statistics excluding v x from the 
selection condition (3.1) which now reads: 

v(t) = v p < v pb = v b (3.7) 
11 



Besides, the variance v b must now exceed the new dynamical border: 

12 

v b > v d = v pd « f p - — (3.8) 

with some empirical factor / p ~ 1 (see Eq.(3.4), and discussion below). 

A typical example of diffusive BF is shown in Fig.2. Both the regular macroscopic 
kinetics of anti-diffusion/diffusion as well as irregular fluctuations around are clearly 
seen. Notice that their size is rapidly decreasing toward the BF maximum. One 
may even get the impression that the motion becomes regular over there, hence the 
term "optimal fluctuational path" [[L4j]. In fact, the motion remains diffusive down 
to the dynamical scale v ~ v <j (3.8). 

Even though a separate BF by itself is fairly regular, the time instance of its 
spontaneous appearance ti, and hence the individual period P, are random in the 
chaotic system. Due to statistical independence of BF under condition (3.6) the 
expected distribution in P is Poissonian (Fig. 3): 

m . ^tmm (3 . 9) 

The principal characteristic of the period statistics, (P), can be estimated as 
follows. From the ergodicity of motion in the iV-dimensional momentum space the 
ratio 

$ = *± = ^4 = %L (3.10) 

This is exact relation (in the limit t run — > oo) where T s is the total sojourn time 
of trajectories within FD (under condition v(t) < t>&) during the whole run time 
t run , and (T s ) is the same per fluctuation. Both ratios are equal to the ratio of N- 
dimensional momentum volume V of the fluctuation at r = to that in equilibrium. 
The ratio $ was also measured during the run. Hence 

(P) = ^ (3.11) 
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The next more difficult step is evaluation of T s = 2T ex from the diffusion equation 
where T ex is the exit (or entrance due to symmetry) time from (or to) FD. A simple 
crude estimate is: T ex ~ Vb/D p = VbC 2 (Section 2). However, the first numerical 
experiments have already revealed that the actual exit time is much shorter, roughly 
by a factor of 1/N 2 . A plausible explanation is in that the most of distribution inside 
FD is concentrated in a relatively narrow layer at the surface of iV-dimensional 
sphere determined by the selection condition v(t) < Vb (3.7). Then, the relative 
width of the layer ~ 1/N implies the observed factor ~ 1/N 2 . Father, the ratio 

®(v b ,N) = v^-^N) (3.12) 

where the geometrical function 

™- 

admits a fairly accurate approximation down to N — 1 (see Fig. 4). 

Collecting all the formulae above, we arrive at our final empirical relation 

(P)~ F ■ 2VbA ° 2 ~ F ■ 2A ° 2 ■ V " N/2 (3 U) 

{P) ~ ¥ ~ F ~W 0(ivy (3 - 14) 

with two fitting factors, A for the layer width, and F for all other approximations 
above. Both factors cannot be united in one because the former enters a new ex- 
pression for the dynamical scale which naturally generalizes Eq.(3.8). Together with 
inequality (3.6) for big fluctuation the new dynamical scale was using for selection of 
purely diffusive BF which are described by Eq.(3.14). The corresponding inequality 
reads (cf. Eq.(3.8)): 

A 12 

v b > v d , v d - — w f p - — (3.15) 

which means that even a small part (A/N 2 < 1) of FD must exceed the dynamical 
scale. 
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Optimization of all empirical parameters was done as follows. The values of two 
factors, B in Eq.(3.6) and f p in (3.15), are not crucial that is the larger they are 
the better for selection of purely diffusive BF. However, this reduces the amount 
of empirical data available. A compromise was found at B — 7 and f p = l which 
leaves 36 runs of 61 done, and 34429 of overall 75053 BF computed with iV = 2 — 10 
for comparison to Eq.(3.14). This was executed in the following way. For each 
selected run with parameters N, C, Vb, and computed (P) and <£> the empirical 
factor F supposed to be a constant was calculated from the first Eq.(3.14). The 
value of parameter A was chosen by minimizing the relative standard deviation 
down to AF/(F) = 0.17. The result, for a given set of data, was A 6. The final 
dependence F(N) is shown in Fig. 5 where the bars are the statistical errors F/y/n 
for each run. 

Coming to analysis of our main theoretical result, the second Eq.(3.14), we first 
remark that it does not describe at all a single trajectory (N — 1). This is because 
we excluded v x b from the selection condition (3.7) (cf. Eq.(3.2)), and thus reduced 
the dimension of phase space to the minimal value, the unity. In this case, a single 
trajectory repeatedly crosses FD with a period P ~ C 2 , the whole diffusion time 
around the phase space torus, independent of FD size. More formally, it follows also 
from Eq.(3.14) since condition (3.6) cannot be satisfied for small iv 

In case of two trajectories (N = 2) the period does not depend on i> 6 , and for 
data in Fig. 5 the ratio (P)/C 2 ps 8.7. Due to fluctuations the actual values of this 
ratio are in the interval 7.4 — 11.0, still not too big for a Big Fluctuation. Apparently, 
this leads to a relatively large scattering of points with N = 2 which also persists 
for iV = 3 too. 

The main, exponential in N, dependence in Eq.(3.14) is readily derived from 
a graphic picture of N statistically independent particles gathering together inside 
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a small domain with probability ~ 1/P ~ v^ 2 Such estimates are known for 
the Poincare recurrences since Boltzmann [RJ. The estimate is especially vivid in 
geometrical picture of iV-dimensional sphere of radius y/v^ considered above. Our 
empirical relation (3.14) considerably improves the simple estimate including a more 
weak power-law dependence which is evident in Fig. 5. 

In our studies described above we fixed the position of BF in the phase space, 
Eq.(2.4). If we lift this restriction, the probability of BF would increase by a factor 
of v b , or by decrease of N by one (N — > N — 1) because now only N — 1 tra- 
jectories remain independent. With the latter change all the above relations would 
presumably still hold true. 

Our main relation (3.14) describes the diffusive kinetics for Vh > Vd, Eq.(3.15), 
when a big fluctuation is not too big. In the opposite case v b «C v d of a very big 
fluctuation, like in Fig.l, the dependence (P(vb)) becomes much simpler 
(see Eqs.(3.11) - (3.13), and (T§): 

< PW > = f * WW) K ^ m (316) 

This is explained by a fast exponential kinetics near BF top (Fig.l) which implies 
the most short exit time T ex 1 and, hence, T s « 2. Indeed, for both BF in Fig.l 
the empirical value of the product = 1.98. 

In conclusion of this Section we show in Fig. 6 the macroscopic kinetics of BF 
entropy, both "exact" (2.6), calculated on the partition of the whole interval 
(0 < p < C) into N p = 401 bins, and that in our approximation (2.7). Both entropies 
were calculated for the same 5 trajectories in one run. A necessary statistics for exact 
entropy was obtained at the expense of a large number n = 4580 of fluctuations in 
the run. For comparison of both entropies we need, first, to adjust the constant So 
in Eq.(2.6). As is easily verified, the Gaussian distribution (2.8) leads exactly to the 



15 



relation (2.7) if the constant 



So = -7; 



2 



In (27re) 



-1.4189 



y/2 



(3.17) 



Approximation (2.7) holds on the most part of BF except a relatively small domain 
near the equilibrium where the distribution in p approaches the homogeneous one. 
The exact entropy (with constant (3.17)) in equilibrium is 



instead of zero in approximation (2.7). The difference is relatively small, the smaller 
the bigger is the fluctuation. In the main part of BF our simple relation of the 
entropy (2.7) reproduces the exact one (2.6) to a surprisingly good accuracy which 
confirms that the distribution in p is indeed very close to the Gaussian one (2.8) as 
expected. 

4 Conclusion: thermodynamic arrow ? 

In the present paper the results of extensive numerical experiments on big entropy 
fluctuations (BF) in a statistical equilibrium (SE) of classical dynamical systems are 
presented, and their peculiarities are analysed and discussed. 

All numerical experiments have been carried out on the basis of a very sim- 
ple model - the Arnold cat map (2.1) on a unit torus - with only two minor, but 
important and helpful, modifications: 

(1) Expansion of the torus in p direction (2.2) which allows for more impressive 
diffusive kinetics of BF out of the equilibrium (Fig.2), and 

(2) Inserting a special (unstable) fixed point for a better demonstration of ex- 
ponential ballistic kinetics as well (Fig.l). Besides, this point was used as a fixed 
position of BF, thus relating our studies of BF to another interesting and important 
problem, the Poincare recurrences (see Eq.(2.2)). 




(3.18) 
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The most important distinction of our approach to the problem was in that we 
have abandoned from the beginning a vague question of the initial conditions, par- 
ticularly a "necessary" restriction of those in statistical physics. Instead, we started 
our numerical experiments at arbitrary initial conditions (most likely corresponding 
to SE), and did observe what the dynamics and statistics of BF were like. In other 
words, we studied the spontaneous BF only. 

What is also important, such a spontaneous rise of BF out of SE as well as 
its subsequent relaxation back to SE can be considered as a statistical macroscopic 
event, even in a few-freedom system like (2.2). The term "macroscopic" refers here 
to average quantities as variance, entropy, mean period, distribution function, and 
the like. 

We consider a particular class of BF which we call the Boltzmann fluctuations. 
They are obviously symmetric with respect to time reversal (see Figs. 1, 2 and 6), 
so that at least in this case there is no physical reason at all for the conception of 
the notorious Hime arrow'. Nevertheless, a related conception, say, thermodynamic 
arrow, pointing in the direction of the average increase of entropy, makes sense in 



spite of the time symmetry [16]. The point is that the BF characteristic relaxation 
time is determined by model's parameter C only, and does not depend on the BF 
itself. On the contrary, the expectation time for a given BF, or the mean period 
between successive fluctuations, rapidly grows with BF size and with the number 
of trajectories (or freedoms), Eq.(3.14). A large ratio of both B = (P)/C 2 1 
is our definition of big fluctuation, Eq.(3.6). A similar result was obtained recently 
|22|j but the authors missed the principal difference between the time arrow and 



in 



thermodynamic arrow. 

A related notion of causality arrow, which by definition points from an indepen- 
dent macroscopic cause to its effect, also makes some physical sense (for discussion 
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sec 



16] and Section 3 above). For Boltzmann's BF considered in the present paper 



the directions of both arrows do coincide independent of the direction of time. The 
latter statement is the most important, in our opinion, philosophical "moral" the 
principally well-known Boltzmann fluctuations do teach us. 

Even though we prefer the discussion and interpretation of our empirical results 
in terms of entropy (S), the most fundamental conception in statistical physics, 
we actually use in our studies another, entropy-like, quantity, the variance v(t) for 
a group of iV trajectories, Eq.(2.5). One reason is technical: the computation of 
v is much simpler while that for S(t) is either very time-consuming in numerical 
experiments (for exact S (2.6)) or approximate (2.7). Besides, for diffusive kinetics, 
we are mainly interested in, the variance is natural variable in which the BF picture 
is most simple and comprehensible. 

Originally, we planned to cover both sides of the BF phenomenon, the regu- 
lar macroscopic kinetics as well as accompanying microscopic fluctuations (noise) 
around. However, our numerical experiments revealed a much more complicated 
structure of the latter as an example in Fig. 2 demonstrates. The dependence v(t) 
looks like a fractal curve on a variety of time scales, from the minimal one ~ 1 itera- 
tion up to ~ C 2 comparable to that of BF itself. This interesting problem certainly 
requires and deserves farther special studies. 

In the present paper the fluctuations in classical mechanics are considered only. 
Generally, the quantum fluctuations would be rather different. However, according 
to the Correspondence Principle, the dynamics and statistics of a quantum system 
in quasiclassics are close to the classical ones on the appropriate time scales of which 
the longest one corresponds just to the diffusive kinetics, and provides the necessary 
transition to the classical limit (for details see ||, |23[). Curiously, the computer 
classical dynamics that is the simulation of a classical dynamical system on digital 
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computer is of a qualitatively similar character. This is because any quantity in 
computer representation is discrete ("overquantized"). As a result the correspon- 
dence between the classical continuous dynamics and its computer representation in 
numerical experiments is generally restricted to certain finite time scales like in the 



quantum mechanics (see two first references [23|). 

Discreteness of the computer phase space leads to another peculiar phenomenon: 
generally, the computer dynamics is irreversible due to the rounding-off operation 
unless the special algorithm is used in numerical experiments. Nevertheless, this 
does not affect the statistical properties of chaotic computer dynamics. Particularly, 
the statistical laws in computer representation remain time-reversible in spite of 
(nondissipative) irreversibility of the underlying dynamics. This simple example 
demonstrates that, contrary to a common belief, the statistical reversibility is a 
more general property than the dynamical one. 
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Figure captions 

Fig.l Mixed kinetics for two big fluctuations of different size. Full/open circles 
show the time dependence of mean variance (v(t — U)) around BF maximum 
at t — tf, upper horizontal straight line is equilibrium while the lower line 
indicates the empirical value of dynamical scale = 0.015, Eq.(3.4) with 
parameter F d rs 1/3. Two oblique straight lines represent the expected fast 
kinetics, Eq.(3.3), and two solid curves do so for the initial diffusive kinetics, 
Eq.(3.5). Run parameters and results are respectively: C — 15, N — 1, v b = 
3.9 x 10- 11 / 6.25 x 10~ 10 (v xb = v pb ), v(0) = 1.96 x 10~ 14 / 3.1 x 10~ 13 , n = 
1971 / 4459, w = 500. Average period between successive fluctuations (P) 

1.4 x 10 7 / 3.5 x 10 6 iterations. 

Fig. 2 Same as in Fig.l for a typical diffusive kinetics (anti-diffusion/diffusion): solid 
curve shows the average over all n = 20259 fluctuations in a run, and wiggle 
line is the same for 28 first fluctuations. Two oblique straight lines repre- 
sent the expected initial diffusive kinetics, Eq.(3.5) with r d = and empirical 
v (em P ) _ q Q45 w hii e the theory (3.15) gives v d = 0.02. Other run parame- 
ters/results: C = 50, iV = 5, v b = 0.0256, w = 10 4 , (P) w 7.7x 10 5 / 8.7x 10 5 , 
and B = 306 / 348; (P)/w w 77 / 87. 

Fig. 3 Histogram of integrated distribution (3.9) for data in Fig.2. Each circle shows 
the number of periods P m > m • AP, m = 0, 1, ... integer; P = n, AP = 

1.5 x 10 5 ; Pmin/w = 1.0027; P max /(P) = 12.63; (P) = 765084. Straight line 
is expected distribution n ■ exp (-P/ (P))- 

Fig. 4 Comparison of directly measured ratio $ ernp , Eq.(3.10), with theoretical ap- 
proximation <& th , Eq.(3.12) for N = 1 — 10: $i = & em p/&th', average over 71 
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runs ($1) = 1.015 ± 0.11 (the standard deviation); the bars show statistical 
errors 1 / ^Jn for each run; the total number of fluctuations in all runs is 127346. 

Fig. 5 Comparison of empirical data for 36 runs, selected from 61 runs computed for 
N = 2 - 10 by the two rules, Eq.(3.6) with B > 7 and Eq.(3.15) with A = 6, 
to theoretical relation (3.14) with the main fitting factor F m , m = 1,...,36 
(see text). Average (F) = 1.51 • (1 ± 0.17) (the standard deviation); the bars 
show statistical errors F m j y/n for each run; the total number of fluctuations 
in 36 runs is 34429. 

Fig. 6 Macroscopic kinetics of BF entropy: lower line is "exact" entropy, Eq.(2.6), 
to be compared with approximation (2.7), middle line; the upper line is the 
same approximation for diffusion theory, Eq.(3.5) with = and empirical 
v (em P ) _ q Q2 j^ un parameters/results: C = 50, N = 5, v b = 0.01, w = 
10 4 , n = 4580, (P) w 3.3 x 10 6 , B = 1314; (P)/w w 329. The number of 
partition bins for calculating (2.6) N p = 401. 



24 



O O f\D 

010 010 



^ - 



T 



o 



7\*) ii i n nm ^ i i m fTThK -N i 

v. \s IFNIK^ I 



I \JLn\Jil 



n.vW/. 



I CB (DgfalGSD I D 1 



Cfr - \—o- 



I-G9H h-OH 



I DG I Cp — I 

00 |- h-OH J O 

I — o— hhoj-i 

I — O — I ll 



o 
o 



cn 



CO 

o 



ro - 



h-OH 



^ - 



o \o\ \-amQ) 



i-CHGO ran 



-H-eei — h 



0> - 



I — O — ^ h-OH 



HOI — I 



00 - 



hfCM 



MO 1 



I O 1 



